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Abstract 

We investigate the effect of the intrinsic spin of a fundamental spinor field 
on the surrounding spacetime geometry. We show that despite the lack of a 
rotating stress-energy source (and despite claims to the contrary) the intrinsic 
spin of a spin-half fermion gives rise to a frame-dragging effect analogous to that 
of orbital angular momentum, even in Einstein-Hilbert gravity where torsion 
is constrained to be zero. This resolves a paradox regarding the counter-force 
needed to restore Newton's third law in the well-known spin-orbit interaction. 
In addition, the frame-dragging effect gives rise to a long-range gravitationally 
mediated spin-spin dipole interaction coupling the internal spins of two sources. 
We argue that despite the weakness of the interaction, the spin-spin interaction 
will dominate over the ordinary inverse square Newtonian interaction in any 
process of sufficiently high energy for quantum field theoretical effects to be 
non- negligible. 

PACS numbers: 04.20.-q, 04.20. Cv, 04.25.-g, 04.25.Nx, 04.40.-b 



1 Introduction 

The Einstein field equations bear the surprising consequence that rotating mass dis- 
tributions tend to drag their local spacetime frame with them as they rotate. This 
well-known effect gives rise to the precession of spinning test masses as well as a 
gravitationally mediated dipole-dipole force (see e.g. [Ij). It has been known for 
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almost a century now that the total angular momentum of a fundamental field, such 
as a Dirac field, is not determined by its orbital angular momentum alone, but only 
the combination of orbital and spin angular momentum. In this paper, we address 
a simple question: Do spinors frame-drag? Or, more specifically, does the internal 
spin of a spinor give rise to frame- dragging effects analogous to that of orbital angular 
momentum? The caveat in the case of internal spin is of course that nothing is actu- 
ally rotating in spacetime, since spin is an abstract, internal property of fundamental 
particles. 

Let us first draw attention to a hole in the current understanding of the coupling 
of intrinsic spin to the gravitational field. For some time it has been known that in 
addition to the standard gravitationally mediated orbit-orbit interaction, there is a 
gravitationally mediated spin-orbit dipole interaction as pointed out by Mashhoon 
[21 El S] 5 Hehl and Ni [5j, Obukhov [B], and others. The existence of such an in- 
teraction follows from a straight forward application of the equivalence principle to 
rotating frames of reference. Consider then the dipole-dipole force between a rotating 
massive object and the intrinsic spin of a Dirac field. As argued from the equivalence 
principle, the force on the spinor will be opposite to the electromagnetic dipole in- 
teraction in that north-pole attracts north-pole and south-pole attracts south-pole, a 
characteristic first noticed by Wald in the context of a spinning body rotating around 
a spinning blackhole [7j. The details of the calculation reveal that the force is medi- 
ated by the gravimagnetic field produced by the rotating mass distribution. On the 
other hand, for Newton's third law to hold in the static case (static so that momentum 
is not carried by the field itself), there must be an equal and opposite force exerted 
by the spinor on the rotating mass distribution. But, what field mediates this inter- 
action? Is it mediated by a gravimagnetic field sourced by the intrinsic spin itself? 
Does intrinsic spin frame-drag? Such a frame- dragging effect has been demonstrated 
for the specific case of a rigid spherical and cylindrical spin-polarized source [H [H [10] , 
however to our knowledge it has not been properly derived from an arbitrary source 
consisting of a Dirac spinor field. If such a source does produce a frame- dragging 
effect, then it leaves open the possibility of a spin-spin dipole interaction. This inter- 
action is distinct from the contact interaction mediated by nonpropagating torsion in 
Einstein-Cartan theory, since it is expected to be a long-distance interactioiJl] coupling 
the intrinsic spins of two spatially isolated spinor fields. 

In this paper we address these questions in the linearized theory of tetrad gravity 
coupled to a spinor field. In particular, in the context of linearized gravity, we show 
that the intrinsic spin angular momentum, or more specifically, the expectation value 
of intrinsic spin, gives rise to a frame-dragging effect exactly analogous to its orbital 
angular momentum. Thus, the total source for the gravimagnetic field is the net 
angular momentum consisting of intrinsic spin plus orbital angular momentum. The 

"'^By long distance here, we simply mean that it is not a contact interaction. General principles 
suggest it should be a dipole-dipole interaction so the force should fall off like r^^. This is considered 
short distance by many comparisons, but it is not a contact interaction. 
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gravitational field at asymptotic infinity must then contain information about both 
the spin and orbital angular momentum, and we show that our recently proposed 
generalization of the the total angular momentum to include intrinsic spin as a surface 
integral at asymptotic infinity [llj does in fact give the right answer in this context. 
Furthermore, under the magneto-static approximation it has the familiar functional 
form as the analogous expression for the ordinary magnetic moment as a surface 
integral following from a straightforward application of Stoke's theorem. 

We then show directly from the gravitational Hamiltonian of an asymptotically 
fiat, static spacetime that, in addition to the orbit-orbit and spin-orbit interaction 
terms of two isolated spinor sources, there is, indeed, a gravitationally mediated spin- 
spin interaction. This is in contrast to what appears to have been the prevailing 
belief that such an interaction could not occur without propagating torsion. In this 
case, the interaction is mediated by the ordinary gravitational tetrad even when the 
torsion is constrained to be identically zero. 

The spin-spin interaction has the same functional form as the analogous electro- 
magnetic dipole-dipole interaction, but with Newton's constant playing the role of 
l/47reo, and with an opposite sign. Since the gravitational interaction is extremely 
weak already, and the dipole-dipole force falls off like as opposed to the fall-off 
of the ordinary Newtonian interaction, it may be argued that in any experimentally 
realizable scenario the spin-spin interaction will be negligibly small as compared to 
the ordinary Newtonian interaction. However, since the spin to mass ratio of fun- 
damental particles is rather large (as compared to, for example, the extremal limit 
of a Kerr black hole), from generic arguments we will show that the spin-spin inter- 
action is expected to dominate at length scales smaller than the geometric mean of 
the Compton wavelengths of the two interacting particles. Thus, as a general rule of 
thumb, we claim that the spin-spin interaction will dominate over the Newto- 
nian interaction in any process of sufficiently high energy for quantum field theoretic 
effects to play a significant role. 

Our conventions are as follows. The Lorentzian metric signature is { — ,+,+,+}. 
Lower case Greek indices {a, /3, fi,!^, ...} represent spacetime indices in the tangent 
space of the four- dimensional spacetime manifold. Lower-case Roman indices near 
the beginning of the alphabet, {a, 6, c, ...}, are spatial indices in the tangent space of 
a three-dimensional spatial hypersurface. Upper-case Roman indices {J, J,K,...} are 
indices in a vector representation of the internal Spin{3, 1) gauge group, and lower 
case Roman indices near the middle of the alphabet, {i,j,k,...}, are spatial indices 
in the fiber. Since we are working in the opposite signature from most quantum 
field theory texts, the Dirac inner product must be defined with an extra i factor so 
that ip = iip'^'^^. There are different conventions in the literature for what the words 
"orbit" and "spin" mean. In this paper, the word spin will always refer to the internal 
spin angular momentum of a spin-| field, and the term orbit will always refer to the 
rotational motion of, say, a rigid rotating object. 
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2 Linearized Gravity in Tetrad Variables 



In order to analyze the geometric consequences of spin in general relativity, it is 
most convenient to consider the linearized theory where many basic concepts such 
as angular momentum and dipole moment that are familiar from flat space theories 
carry over to curved spacetime. In order to couple spinors to the gravitational field, it 
is necessary to resort to the tetrad formulation of gravity where the dynamical metric, 
gfj^i, is replaced by the dynamical tetrad e^^ related to the metric by g^u = rjije^ ^ 6"^^, 
where rju = diag{—l, 1, 1, 1) is the canonical bilinear form. In the linearized theory, 
the tetrad is expanded about a fiducial flat metric °e^^, which for convenience we will 
set to be the canonical tetrad ^e^ ^ = 6^ ^ so that the tetrad and its inverse become 

e'^ = 6'^ + e', e^j = 5^j-e-j. (l) 

The perturbations e are assumed to be small enough that we consider only expressions 
that are first order in e. By convention, in the linearized theory, we will raise (lower) 
indices using the fiat metrics tj^'^^tikl) and r]"^{r]^u) and convert between internal and 
spacetime indices with fj_{S''j). Ordinarily the index placement of is immaterial 
since the object is itself the map between the tangent space and the internal SO {3, 1) 
vector space - however, since the fiducial fiat tetrad now plays that role, the index 
placement of does matter. Specifically, the relation between e and the ordinary 
metric perturbations is h^^, = 2 £(^y) where g^^ = rj^^ + h^^ and following our conven- 
tions, e^y = Tjjjd^ ^s'^ y. We will adopt the convention that antisymmetric component 
fields will be denoted with an underbar, and symmetric components will be denoted 
with an undertilde so that e\ = ej ^ ~\ "^i^^ ej ^ = and ^ = ri^°'e(^afj.)- 

Since we wish to emphasize that a spin-spin interaction arises even without a 
propagating torsional field, we will work in the second order formulation of Einstein- 
Hilbert gravity formulated in terms of tetrads. For more on the low-energy effects of 
torsional theories, see for example [13] . Given a basis of coframes, the torsion-free 
spin-connection coefficients can be written 

T''^[e] = 2e[^l^le-^][,,^] + e^'^e'^CM^.e^' [a,/^] (2) 

which in the linearized theory reduces to 

r^-^^ = 2r/^[^£^][,,^] + v^KV^'v^'-'e'^M • (3) 

Writing the curvature R^"^ = 2(9[^r^"^,^] + xifJ'^'^u]) and the Ricci tensor R^ ^ = 
e^jR^-^^u, the Einstein equations are 

G', = R'^-h'^R = 87rGr\. (4) 

As usual, by choosing a gauge in which the divergence of the perturbation of the 
tetrad vanishes, e^o-''^ = 0, the linearized Einstein equations take the form 

-De'^ = 16nGT'^. (5) 

where = ~ k^fi^ e = 6j ^. 
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3 Fermionic Matter Sources 

We will now take the source to be a Dirac field. The total action for the system is 

S = Sg + S^ (6) 



IGnG _ 



(7) 



Here we have included a coupling constant a to account for different conventions. 
Normally, the constant can be taken to be a = —1. 

Variation with respect to the spinor field gives the Dirac equation 

e'^n^Df^ip = mip (9) 

where D^ip = d^ip + |7[/7j] ^^'^ ^l i'- Variation with respect to the tetrad yields the 
Einstein-Cartan equations: 

R'^ - ]^e\R = 8nGr'^ (10) 

where the stress-energy is given by 

T^^ = -a e-^" ^ {i^l^aD^)^ - Df^^^^a^ip) . (11) 

For future reference, it will be useful to mention the following identity which follows 
from the extra six components of the fixed points of the variation of the action with 
respect to e^^, but also follows from the Dirac equation: 

I (^Tb^MlV' - D^^^la]^) = -'-ea^DaJj,'^ (12) 



(5) ^ 

where Jp = ip'y^'ypip is the axial current. 



4 Multipole Expansion of the vector potential 

In this section we review the generic procedure for calculating the vector potential 
from an current distribution in the monopole-dipole approximation (for the electro- 
magnetic case, see e.g. [H]). The generic wave equation 

- DA^ = kJ^" (13) 

can be inverted using retarded Green's functions to yield 

Af^ix) = — I d'x'j^^^^^ (14) 
^ ^ Arc L r - r' 
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where the integral is over the past hght cone. Under the assumption that the current 
distribution is locahzed near the origin, one can expand the integrand in powers of r 
so that \r — r'\~^ = 1/r + f ■ r'/r"^, yielding 

A''{x) = ^(- [ dVj^(x') + ^ / d^x'x"'J^ix')£x'] , (15) 



where = x^—6q x° is the spatial radial vector and f'^ = r'^/r. Under the assumption 
that the current is localized and conserved, d^J^ = 0, one can derive the following 
identities: 

f d^x'J" = [ d^x'J°+ [ d^x'x'^'dtJ^ (16) 
(x'^V')) = / t/V^JV^^JV / t/Vx'V^9tJ°. (17) 



To simplify the expressions, we make the standard "magnetostatic" assumption that 
the charge distribution is conserved in time so that we have the additional constraint 

dtJ^ = daJ^ = 0. (18) 

Under this assumption, the field distribution is static so the integral over the past light 
cone can be replaced by an integral over a constant time slice. Then, by splitting the 
second integral on the right-hand side of f[T^ into its symmetric and antisymmetric 
components, the expression reduces to 



where 

Q = j d^x'J° (20) 

is the total charge, 

p = [ d^x' r'.f (21) 



is the electric dipole momenio, and 

M^" = I d^x'x'^^'r^ (22) 



E 



is the relativistic generalization of the magnetic moment whose spatial components 



are related to the ordinary magnetic moment by (m)" = ^e°'i,cM^'^. Splitting this into 



^As we have written expression ()19p . half of the electric dipole momentum appears in the rela- 
tivistic angular momentum, explaining the explicit ^ in the expression. 
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spatial and time components, we have 



^° = ^ ? + ^ (23) 



A' = f (24) 

Focusing on the spatial components of the vector potential alone, we rewrite the above 
in vector notation as 

A- ^ ^ X ^ (25) 
which gives rise to the dipole magnetic field 



B = V X A = — 

4:71 



3f(fm) — m Stt ^o, , 

1^13 + 



(26) 



The delta-distribution term in the magnetic field is the well-known Fermi contact term 
that gives rise to the hyperfine interaction in the case of ordinary electromagnetism. 
For future reference, it is required so that a surface integral expression for the dipole 
moment agrees with a bulk integral expression. Specifically, given the form of the 
vector potential A, it is easy to see that 



X A B?dn = m . (27) 



Using Stoke's theorem to yield J^^^ 'V X A d^x — R"^ J f X A dQ, we see that the 
bulk and the surface integrals only agree if the delta-distribution is included in the 
magnetic field. 



4.1 The magnetic dipole field of a spinor 

To derive the magnetic dipole field of a spinor, we take the vector potential Af^ to be 
the ordinary electromagnetic vector potential so that k, = fiQ, and we take the current 
to be the electromagnetic current of a massive Dirac field with charge q: 

J'^iqfn.rac-iQh"^- (28) 

We first notice that subject to the Dirac equation, j^Dnip = mip where = — 
iqAfj,, the Dirac current splits into two, separately conserved currents: 

j Dirac ~ JkG + Jspin (29) 

= ^ {^D^i^ - {D^i^)i^) + (V;[7^ tIV') ■ (30) 

This splitting is often referred to as the Gordon decomposition of the current. We em- 
phasize this step because a similar decomposition will allow us to separate the orbital 
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and spin components of the corresponding gravitational current, the stress-energy. 
Using this identity, upon integration by parts the relativistic magnetic moment be- 
comes 



M"'^ = q 



(31) 



Note that the term coming from the spin current contains no derivatives, so it de- 
pends only on intrinsic spin properties of the Dirac field. As usual, to examine the 
nonrelativistic limit of the Dirac equation, we work in a representation where the 
Dirac matrices take the form: 



7 





' 1 " 


7* = —i 


" a^' 




"01" 


i 


-1 





75 = 


1 



(32) 



In the nonrelativistic limit, the Dirac field splits into large and small components: 



(33) 



and the large components dominate, so that, for example, the relevant currents can 
be approximated by 

iip'^^ip ^'^(^ —2e''\^'^(j^^ z?/'757Y ~ — . (34) 

In this limit, the magnetic moment becomes 



2 2m 



d^x —e"^^ - (x6D,$"f)$) + 2x [ d^x' i$V"<l> 

2 ./v 2 



(35) 

To understand this better, we compare this result to the Noether charge associated 
with angular momentum given by (see e.g. |llj ) 



Q{ij} 



(36) 



The nonrelativitic limit of this expression yields the angular momentum expression: 



'tot) 



d?x Y^^bc {'^^x''D^<^ - (a;^D"$^)<l>) + J d^x ^^V"*!' . (37) 

The first integral is clearly the expectation value of the orbital angular momentum, 
and the second is the expectation value of the spin. Comparing these expressions, we 
obtain the expected form for the magnetic moment: 



rriEM = -^{L + 2S) . 
2m 



(38) 
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It is therefore clear from the discussion leading up to fl26|) that the internal spin of a 
localized Dirac field produces its own dipole magnetic field characterized by a dipole 
moment of twice the spin angular momentum. 

For the gravitational interaction, the procedure is very much analogous to the 
above derivation. In this case, the linearized perturbations to the tetrad will play 
the role of the gauge field, which is sourced by the stress-energy as the conserved 
current. A corresponding "gravitational Gordon decomposition" of the stress-energy 
will allow a separation of orbital angular momentum and spin. In place of the Dirac 
bilinear ^'y^^^'^'^^ip as the intrinsic spin source for the magnetic field, the axial current, 
will play the role of the intrinsic spin source for the gravimagnetic field. 



5 The dipole expansion of the gravimagnetic field 

The perturbation to the tetrad in the linearized theory satisfies the same set of equa- 
tions as the vector potential given above. That is, we can think of set of four 
vector potentials labeled by the internal index J. The coupling constant is k = SttG. 
Thus we have the immediate solution: 



or separating space and time components. 



= -2G [9L + tjf] (40) 



- 2g '^ (41) 



where by analogy we have defined 

Ml, = jj'x'xlr',]. (42) 

Let us now focus on the time components . Analogous to ordinary gravito- 
magnetism, in a low energy, low velocity limit we expect these terms to dominate. 
Since the resulting field equations in this gauge are very similar to the electromagnetic 
field equations in the Lorentz gauge, as in ordinary gravitomagnetism we define the 
gravito-electromagnetic vector potential An = . In the nonrelativistic limit, the 
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charge reduces to the expectation value of the energy: 

g° ^ -a{E) = -a [ d^x'- {<!>Udt<!> - (43) 

and the dipole p° is the dipole associated with the energy density. 

We now focus on the magnetic moment term M^^,. As in the case of electromag- 
netism, the first step is to rewrite the matter current so that the internal spin of the 
particle is more transparent in the current. To do this, we employ the identity fll2p . 
which allows us to write the stress tensor in a more convenient form analogous to 
the Gordon decomposition of the Dirac current (for a more in-depth review of the 
Gordon decomposition of the gravitational stress-energy and spin-currents, see |15j): 

r', = -a Q i^Pi'D,^ - D.i^y^) + Ig^/'^D^ ) (44) 

In the linearized limit, the components M°j^ in the magnetic dipole contribution of 
the multipole expansion reduce to 

= -a I d^x'^ (V^7°a;[^9,]7/; - (x[^9,]V^)7V) - ^eo^.^V^TsT^^ • (45) 

We recognize this as proportional to the conserved Noether charge associated with 
angular momentum: 



M' = -iQi,u} . (46) 



Thus, from our previous calculation, we see that the spatial components of the grav 

1 

2 



imagnetic moment, (^n^)" = ^e°'hcM^'^, are precisely one-half of the total angular 



momentum: 

m, = -|(L + 5). (47) 
The resulting gravimagnetic vector potential has the form of a dipole potential, 

A=2G'I:^. (48) 



which gives rise to the dipole gravimagnetic field 



B = V X A = 2G 



3f(f • m„) - m„ Sir ^n, , 



(49) 



In the special case where the particle is localized and at rest, so that the angular 
momentum vanishes, the spin of the particle itself produces the gravimagnetic dipole 
field: 

''-fi^4^.'-^SSH.l (50, 



B'P''' = -Ga 
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6 Computing the internal spin angular momentum 
in the linearized limit 



We recall that in magnetostatics, the magnetic moment of a system can be computed 
by a surface integral at infinity as follows: 



m = lim ^ f r'^dVt f x A . (51) 
'52 



2k 



Since the monopole contribution is zero, and the higher order terms in the multipole 
expansion fall off faster than r~^, the integral picks up precisely the dipole contribu- 
tion of the magnetic field. In a recent paper, we have constructed a similar expression 
for the total (spin + orbital) angular momentum of an asymptotically fiat spacetime 
in the tetrad framework, and related it to the Noether charge corresponding to an- 
gular momentum in the axisymmetric case. We will now show in the linearized limit, 
that this integral gives the analogue of (15T|) for the gravimagnetic field. 

The expression for the total angular momentum derived in pLlj is the following 
(here we have normalized the expression in [TT] by an overall minus sign): 

4i = L^' + S^J (52) 
L'^ = ^ [ euKLtKiiJiie' Ae-')Aco'''^ (53) 



s 



S'' = e^^^.Ajf A . (54) 



E 



In the above expression, k = SnG, and the vectors K^^"^^ form a six-dimensional basis 
for the boost and rotation Killing vectors of the fiducial flat metric "(7^1^; the internal 

boost and rotation parameters XY/^ are related to Killing vectors by 

C^^uy'e' = -i\^''^yK'e^. (55) 

In the Cartesian basis where "^e^^ = 6^ ^, the Killing basis vectors can be written 
K^~j = —X'^'^j Xu and the internal boost and rotation parameters are given by 

Xiji-. = Sj^S^ -^/f 71 = ^i^i ~ 5\5i. In addition to this, one also has the linear momen- 
tum given by (again correcting by an overall minus sign) 



1 

4fc 



ewL^;f{/}(e^ Ae-^) Aw^^ (56) 



where K^^^ is a translational Killing vector for the fiducial fiat metric given in the 
Cartesian basis by K^^^ = SItt-- 

The key property of the spin and orbital angular momentum generators is that 
they satisfy the algebra of the spin enlarged Poincare group under the Poisson bracket. 
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The spin enlarged Poincare algebra is the subgroup, (5(°e) C Spin{3, 1) xi Diff4^, of 
the local gauge group of Einstein-Cartan gravity that preserves the fiducial flat tetrad. 
The group has the structure of Spin{3, 1) (g) (5*0(3, 1) x M^'^) together with a two-to- 
one mapping between the Spin{3, 1) and S0{3, 1) subgroups that ensures that the 



group preserves °e. 



We will now evaluate this expression for the spatial angular momentum compo- 
nents JjQ^ = \^^jkJtof We will work under the assumption that the Killing vector 
is rotational, and the gravitational field is static. In this case, the orbital angular 
momentum integral yields 

= 11 r-K^^' (£o,,« - lo., J r'dn (57) 

= -llj">^iH^£o'-^o'')r'dn (58) 
and the spin contribution is 

S'^ = -\j f'^X^ (£o^ +fo'') ^^dVL (59) 
Adding these together we obtain 

Jlot = l^d^ + S^) = -^ [ r'^^yX^^^^eo'-r'dn (60) 

(f X Af (61) 



SttG 
3 



52 

Bf . (62) 



8nG 

From the form of the gravimagnetic field in the nonrelativistic limit (H9|l . 

4 = -«(X + 5)\ (63) 

Using the standard convention a = — 1, this gives precisely the total (spin+orbital) 
angular momentum of the static current source. 



7 The Spin-Spin dipole interaction 

We can now partially answer one of the paradoxes discussed in the introduction. The 
well-known spin-orbit interaction as lauded by Mashhoon, Hehl and Ni, and others, 
gives a dipole-dipole interaction between the orbital angular momentum of a rotating 
matter distribution and the intrinsic spin of an elementary particle such as an electron 
or a neutrino. Opposite to the electromagnetic interaction, like poles attract. The 



12 



force that the rotating mass distribution exerts on the intrinsic spin of the spinor is 
mediated by the well-understood frame-dragging effect of orbital angular momentum. 
On the other hand, if Newton's third law is to hold in the static Newtonian limit, the 
spinor must also exert an equal and opposite force on the rotating mass distribution. 
But what is the field that mediates this force? From the above discussion, the answer 
is now clear - the intrinsic spin itself gives rise to a frame- dragging effect, which 
gives rise to the force that the spinor exerts on the rotating matter distribution. On 
the other hand, this opens up the possibility of an additional spin-spin interaction in 
addition to the well-known spin-orbit interaction. Furthermore, since the spin to mass 
ratio of fundamental particles is abnormally large, we expect the spin-spin interaction 
to be the dominant gravitational interaction at short range. We will now derive the 
weak field Hamiltonian and show that there is a characteristic spin-spin interaction 
term. 

7.1 The weak field Hamiltonian 

We wish to derive the interaction Hamiltonian in a gauge invariant manner. Gener- 
ically, due to diffeomorphism invariance, the total Hamiltonian in a spatially closed 
universe is identically zero. However, in the asymptotically flat context the Hamil- 
tonian is not zero, and it is the additional boundary terms which give rise to the 
nonzero Hamiltonian. Under the assumption of a timelike Killing vector, appropriate 
in the static case, the boundary integral can be related to a bulk integral of the matter 
content, which in our case will give rise to the interaction Hamiltonian. In previous 
sections we have made the magnetostatic assumption which at the linearized level 
implies the existence of a timelike Killing vector. We will work with the covariant 
phase space consisting of solutions to the full set of equations of motion. As shown 
in [TT] the relation between the Hamiltonian as a boundary integral, and the bulk 
integral under the assumption of a timelike Killing vector is given by 



tot 




(64) 




a 



6 



(65) 



(66) 



where the induced metric volume on the three-surface given by a^ ' = 

e"^ A A e^). Employing the on-shell Dirac equation this becomes 



iti^euKie^A 



(67) 




(68) 
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and we describe each term separately: 



• Hq-. This is the free field Hamiltonian density with an overall modulating factor 
to account for the gravitational redshift giving rise to the Newtonian potential: 



H'o^aJ d^x (1 + e) (^^ 



{ijj^°'daijj — daip^'^ip) — mi/jip 



(69) 



• H^: This term gives rise (in our choice of gauge) to the couphng of the orbital 
angular momentum of the particle to the external gravimagnetic field. In other 
words, it gives rise to the Bg^t " L coupling, and it is given explicitly by 



(70) 



H}^. This term couples the external gravimagnetic field to the internal spin of 
the particle and gives rise (in our gauge) to the B^xt • S coupling. It is 



-a 



j^d?x ^^{7°,ro}v^. 



(71) 



ff(3) -|- ff(4) - These terms depend explicitly (in our gauge) only on and 
therefore couple the pressure and shear of an external source to the spinor. In 
the gravi-magnetic approximation it is assumed that the energy-momentum and 
angular momentum will dominate over the pressure and shear terms. Thus we 
will neglect these terms. For completeness, they are given by 



^(3) + ^(4) 



-a 



(fx 



2 



(72) 



To determine the precise form of the spin-orbit and spin-spin coupling, we make 
the following assumptions. Suppose we have two spinor sources ipi and 1^2 highly 
localized in regions Ui and U2, respectively, that are widely separated by an average 
distance r. The interaction Hamiltonian we are interested in comes from the terms 
if*^^) -I- i?^^). We will assume that the hnearized gravitational field coupled to V'l 
is sourced entirely by ip2 and vice versa. Thus, we will ignore the self-interaction 
of the field. The symmetric components e^^ arc determined in the linearized limit 
(in the Lorenz gauge) by the fermionic sources as shown above. The antisymmetric 
components are fixed by a choice for the internal Spin{3, 1) gauge. The calculations 
simplify if we choose a gauge such that £5^ = 0. With this gauge choice the tetrad 
perturbations takes the form 









^00 




= 2G 
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(73) 
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where it is understood in this expression that f is the unit vector pointing from the 
source to the target and r is the distance between the two as measured using the 
fiducial flat metric. 

Focus now on the first term in the Hamiltonian: 

ifW = -4G'a/" rf^xi ma- ( X - (^/'7°V^/'- VV-tV) I +(1^2) (74) 

where r2i = r\ — r2- Under the assumption that the regions Ui and U2 are small 
compared to the distance separating them, we can make the expansion in the region 

f.x _.r.-r._/.{f.-rr)_ ^^^^ 



r2i| r2| 
Using this approximation, and assuming the expectation value of the linear momen- 
tum of both particles is zero (which is necessary in the static case), after some ma- 
nipulation the expression for H^^^ reduces to 

- -2Ga i^irn,.f,)iL,.r,)-m,.L,) ^ ^ ^^^^ 

We now turn to the second term in the Hamiltonian, H^'^\ From the identity 
{7"^, 7['^7'^1} = 2ie^"^'^^757/,, and the linearized form of the spin connection, F/jo- = 
SjjEjj-^^u — €jj^^, the expression reduces to 

H^'^ = -a I V*^'=V^i757,^iF,,o + (lf^2) (77) 

^ ci^xi e^^'=9,£o^ (^Y^i757fcV^i) + (1^2) (78) 

^ -aB2-Si + (1^2) (79) 
where B2 is the gravimagnetic field produced by ip2, given by ( 149|) . Thus, in total we 

ff(^) ^ -2Ga +(1^2). (80) 

Putting both these terms together and inserting the form of the gravimagnetic 
moment m = —^{L + S), we have the total gravimagnetic interaction 

Hii)^H^2) ^ ,^^2 3 {{L, + 5i) ■ r) {{L^ + S^) ■ f) - {L^ + S^) ■ (L^ + 5^) 

2 3 ( J*°* ■ f) ( J*°* ■ f) - Jf* • J*"* 



2Ga' 



In addition to the orbit-orbit and the spin-orbit terms, we have a new spin-spin 
interaction given by 

n - orr.^ HSi-f){S2-r)-Sr-S2 

-'^ spin— spin ^wCt ^ . \^^) 
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As expected, this interaction has the usual form of a dipole-dipole interaction, but in 
contrast to the electromagnetic interaction, north-pole attracts north-pole, and south- 
pole attracts south-pole. A similar interaction appears in the context of propagating 
torsion theories fi6\ [13], where the torsion mediates the interaction, as well as in 
the context of Maxwellian gravity [17]. Perhaps because of the former result, the 
prevailing belief in the community appears to be that such a long-range interaction 
could not occur without propagating torsion, or could only be manifest as a contact 
interaction with nonpropagating torsion. In our case, the interaction comes from an 
ordinary gravitomagnetic effect derived from the coupling of tetrad gravity to a spinor 
field, and it is not mediated by torsion. 

8 When is the spin-spin interaction non-negligible? 

The Newtonian gravitational interaction between two fundamental particles is exceed- 
ingly small at large distances, and since the spin-spin force falls off like ^, this force 
will be even smaller. However, as the particles come closer together eventually the 
magnitude of the spin-spin interaction will dominate over the ordinary attractive ^ 
interaction. Fundamental particles have the peculiar property that their intrinsic spin 
to mass ratio is rather large, as compared to, for example, the extremal limit on the 
angular momentum of a Kerr black hole. Thus, the spin-spin interaction may domi- 
nate at distance scales much larger than one might ordinarily expect. We can make 
an heuristic argument to determine the distance scale at which the force becomes 
non-negligible. This will be true when the Hamiltonian for the ordinary Newtonian 
attractive force is of the same order of magnitude or less than the spin-spin inter- 
action Hamiltonian. Using the minimum values of the interaction^, the spin-spin 
interaction will dominate over the Newtonian interaction within the range given by 





< 




r 







or, including the factor of c. 



where Ac = and Ac = are the (reduced) Compton wavelengths of the inter- 
acting particles. Using the explicit form for the tetrad perturbations, [731 it is clear 

■^Of course, in contrast to the r^^ Newtonian force law, the r^^ spin-spin force can be both 
positive and negative, and therefore has the potential to cancel itself out at the microscopic level. 
By "dominate" interaction here, we simply mean the interaction which we can never justifiably 
neglect without additional assumptions. 



rriic 

For spinors the maximum magnitude of the spin is s = |, so we obtain 
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that even at these length scales we have \s^y\ ^ 1, so the linearized approximation 
is still valid near the Compton scale. Thus the spin-spin interaction is expected to 
be the dominant interaction between fundamental particles at length scales smaller 
than the geometric mean of the Compton wavelengths of the interacting particles. 

The spin-spin interaction involves the internal angular momentum, which also 
couples to the electromagnetic field to yield a dipole-dipole coupling. Since the elec- 
tromagnetic interaction is much stronger than the gravitational interaction, the effect 
will be most prominent in the dynamics of weakly interacting neutral particlefl Tak- 
ing the upper limit on the lightest neutrino to be around 1 eV , the scale at which the 
spin-spin interaction dominates is around 10~^ cm. At these length scales, we are well 
within the bounds of validity of the linearized theory. On the other hand, the Comp- 
ton wavelength effectively sets the scale at which quantum field theory effects become 
non-negligible since the Compton wavelength is the wavelength a photon would need 
to create a particle of mass m out of the vacuum. Evidently, the regime where the 
spin-spin interaction becomes significant coincides with the regime where quantum 
field theoretical effects become appreciable. Thus, we propose the following "rule of 
thumb" : 

In any process of sufficiently high energy for quantum field theoretical 
effects to be appreciable in the interaction of two spin-^ particles, the 
gravitationally modulated spin-spin interaction is expected to be of equal or 
greater significance than the ordinary attractive Newtonian gravitational 
interaction. 

Alternatively, one could say that quantum field theoretical effects are likely to be the 
most prominent in the dominant interaction in the regime where quantum processes 
are non-negligible. Following this line of reasoning one may conclude that quan- 
tum effects should be the most prominent in the spin-spin interaction. To support 
this argument we can compare the situation with that of quantum electrodynamics: 
In classical electrodynamics, the scale at which the spin- spin interaction potential 
/io(Aii ■ ^^2 — 3^1 ■ r fj,2 ■ r)/4:iir^ dominates over the charge potential gig2/4vreor is 



not surprising that the most prominent quantum effects of QED are manifest in the 
spin-spin interaction, namely in the correction of the gyromagnetic ratio. Although 
quantum gravity effects are expected to be dramatically smaller, following this same 
line of reasoning one might expect that some of the most prominent effects of pertur- 
bative quantum gravity or semiclassical gravity with backreaction considerations will 
involve the spin-spin interaction. 

Strongly interacting particles are likely to be poor candidates for phenomenology of this interac- 
tion because they are typically composite particles made of electrically charged quarks. For example, 
despite being electrically neutral, the neutron has a small ordinary magnetic moment coming from 
its composite nature. Despite being small, it is likely that the electromagnetic spin-spin interaction 
between neutrons would dominate over the gravitational spin-spin interaction at all length scales. 



also classically given by the geometric mean 




In light of this, it is 
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As a final comment, we note that the primary purpose of this paper is to demon- 
strate the existence of a generic spin-spin interaction. It remains to be seen if this 
interaction has any apphcation to the standard model of particle theory, cosmology, 
or astrophysics. To answer these questions, a more realistic model is necessary. Such 
a model would necessarily take into account the effects of the weak and strong inter- 
action, and since these interactions are generally much stronger than the gravitational 
interaction, and sine each has its own spin-spin interactions associated with it, it may 
seem that weak and strong effects will override the spin-spin gravitational interaction 
in the regime where it begins to dominate over the Newtonian potential. However, a 
simple heuristic argument shows that there may be a window where the gravitational 
spin-spin interaction may be the dominant, albeit small, interaction. As a simple 
model consider the wave equation for a massive vector boson (which can be used to 
model the gauge bosons of the weak force in a linearized approximation) : 

(□ - {^Jic/Kf)A^ = -i^r (86) 

where ii is the mass of the gauge field. As shown previously, using a retarded Green's 
function this relation can be inverted to give the Yukawa vector potential 

A'^^x) = / ^ ^-i^ (87) 

Ja_ \r -r' \ ^ ^ 

where = ^ is the Compton wavelength of the massive boson. From this it is clear 
that the weak effects of a massive gauge mediated force fall off rapidly outside of the 
scale set by the Compton wavelength, r ~ A^. For weakly interacting particles, such as 
neutrinos, the gauge bosons have a mass of around 100 GeV corresponding to a length 
scale around 10~^^m. As we have seen, the spin-spin interaction for particles with 
masses around leV dominates at around 10^^ m. So there is a large window where 
the spin-spin interaction could potentially dominate over the Newtonian interaction 
and not be washed out by electroweak phenomena. It should be stressed that this is 
a heuristic model only, and more accurate predictions would require a neutrino mass 
extension of the standard model together with semiclassical gravitational couplings. 
More thorough investigations of this sort are currently in progress. 

9 Concluding Remarks 

In this paper we have investigated the gravitational effects of a Dirac field as a source 
of matter, focusing on the effect of the internal spin-angular momentum on the space- 
time geometry. We have shown that in the context of the linearized theory, despite 
the absence of a physical rotating mass distribution, the intrinsic spin of an Dirac field 
gives rise to a frame-dragging effect precisely analogous to that of ordinary orbital an- 
gular momentum. This resolves some of the paradoxes mentioned in the introduction. 
In particular, in the Newtonian limit, Newton's third law is regained for the spin-orbit 
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interaction because the orbiting mass experiences a force mediated by the gravimag- 
netic field produced by the spinor itself. This yields the surprising consequence that 
with two spinor sources, the gravimagnetic field mediates a long-range spin-spin in- 
teraction analogous to the dipole-dipole interaction of ordinary electromagnetism. It 
remains to be seen if such an interaction has any astrophysical or cosmological con- 
sequences. The effect will be most pronounced in the interaction of a neutral particle 
with another spin source since otherwise the corresponding electromagnetic interac- 
tion will dominate over the gravitational interaction. Thus, we expect that such an 
interaction may be significant in the dynamics of neutrinos as they interact with other 
spinor fields. The consequences of such interactions are currently under investigation. 
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